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Abstract 

We recalculate the two-loop beta functions in the two-dimensional Sine- 
Gordon model in a two-parameter expansion around the asymptotically free 
point. Our results agree with those of Amit et al., J. Phys. A13 (1980) 585. 



In this paper we recalculate the two-loop beta-function coefficients in the two- 
dimensional Sine-Gordon (SG) model in a two-parameter perturbative expansion 
around the asymptotically free (AF) point. The study of the SG model in the vicinity 
of this point is especially important since this region is used in the description of the 
Kosterlitz-Thouless (KT) phase transition in the two-dimensional 0(2) nonlinear o- 
model, better known as the XY modelQ. This was the motivation of the authors of 0, 
who have undertaken a systematic study of perturbation theory in a two-parameter 
expansion around the AF point. They calculated the renormalization group beta- 
functions up to the two-loop coefficients. The beta-function coefficients were also 
calculated in ]3| by a completely different technique based on string theory. The 
results found in || differ from those of at the two-loop level. The question of two- 
loop beta-function coefficients were considered also in M for a class of generalized 
Sine-Gordon models. The results, when specialised to the case of the ordinary SG 
model, agree with those of M, but disagree with those of [0. In |5|] the short 
distance expansion of some Sine-Gordon correlation functions were calculated using 
conformal perturbation theory. This allowed the extraction of the one- and two- 
loop beta-function coefficients around the AF point. The resulting two-loop beta- 
functions differ from all the previous results. 

In view of the role the Sine-Gordon model is playing in the description of the 
KT phase transition it is very important to know the correct two-loop beta-function 
coefficients. The purpose of the present paper is to show that, in fact, the two-loop 
results of Amit et al. are the correct ones. We show this first by comparing the SG 
beta-function to known results in the chiral Gross- Neveu model |3| , which is known 
to be equivalent to the SG model at its AF point. We also check the beta-function 
coefficients by considering the renormalization of 2n-point functions of exponentials 
of the SG field. 

Following m we consider the Euklidean Lagrangian 

C = \d^d^ + + J|L [1 - cos(A,0)] , (1) 

where mo is an IR regulator mass and a is the U V cutoff (of dimension length) . UV 
regularized correlation functions are calculated by using 

G Q (x) = ^K (mov^T^) (2) 

where K is the modified Bessel function, as the propagator. Our strategy is 
slightly different from 0, who really considered the renormalization of the massive 
SG model ([I]) of mass itiq. We treat mo as an IR regulator mass and consider IR 

1 For a review of the Sine-Gordon description of the Kosterlitz-Thouless theory, see Q 
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stable physical quantities for which we can take the limit mo — > already at the UV 
regularized level (before UV renormalization). 

The model ([I]) is renormalizable in a two-parameter perturbative expansion 
around the point corresponding to the couplings a = 0, j3 2 = 8n. Writing 

(3 2 = 8n(l + 6 ) (3) 

the two bare expansion parameters are a and <5 and physical quantities can be 
made UV finite by the renormalizations 

a = Z a a; Z Q = l+ gi 5£ + a 2 (g 2 £ 2 + g 2 £) + 5 2 (g 3 £ 2 + g 3 £) + ..., (4) 
l + 5 = Z^(l + 5); Z 4> = l + f 1 a 2 £ + a 2 5(f 2 £ 2 + f 2 £) + ..., (5) 

where a and 5 are the renormalized couplings and £ = In fia with /i an arbitrary 
renormalization point. The dots stand for terms higher order in perturbation the- 
ory and the numerical coefficients g±, f\ etc. can be calculated by renormalizing 
correlation functions. The results of Amit et al. are 

fi = 9i = -2, / 2 = -Jj, 92 = 93 = 0, (6) 

those of || and || are 

fi — 32) 9i — ~ 2, /2 = —33, 92 = —35, #3 = (7) 
and finally f5[ found 

fi = §5, 9i = -2, /a = -gg, £2 = -53, #3 = 0. (8) 

We see that the one-loop coefficients are the same but not all two-loop coefficients 
agree. The subject of this paper is to recalculate these numbers. 

The renormalization group (RG) beta-functions can be calculated by solving the 
equations 

Va = V5 = 0, (9) 

where, as usual, the RG operator is defined by 

d d d 

V = ~ a da + ^ a °' So ^da~ + P 5 ( a °> S ^dS~' ( 10 ' ) 



One finds 



p a = -g iao S - g 2 a 3 - g 3 a 5l + . . . , (11) 
Ps = fi<4 + (/1 + fiHSo + . . . (12) 
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It is well-known that, in the case of several couplings, the higher beta-function 
coefficients are not all scheme independent. Indeed, considering the perturbative 
redefinitions 

a = a + cia 5 + . . . 5 = 5 Q + c 2 al + . . . (13) 

one finds that in addition to the one loop coefficients f\ and g\ only the following 
two two-loop coefficient combinations are invariant. 

g 3 , J = 2g 2 - f 2 . (14) 

The RG analysis with two couplings can be made similar to the case of a single 
coupling by changing the variables from «o and So to the pair Q and So, where Q is 
a RG invariant (solution of the T>Q = equation), given in perturbation theory by 

Q = ha 2 + 9l 5 2 + 2g 2 a 2 5 + F 2 5 3 + . . . , (15) 

where F 2 = |g 3 — |gx + ^ J. Now Q, being a RG invariant, can almost be treated 
as if it were a numerical constant and So as the 'true' coupling. The beta-function 
in these variables is 

(3{5 , Q)=Q + 25% + A Q$o + B5I + ..., (16) 

where A = 1 - J/ f x and B = 2{A - g 3 )/3. 

It is well-known that the SG model can also be formulated in terms of two 
fermion fields, interacting with a chirally symmetric current-current interaction [|TJ. 
A special case of the two-fermion model corresponds to the S*?7(2)-symmetric chiral 
Gross-Neveu model. This correspondence is evident in the bootstrap aproach, since 
the SG S-matrix in the limit (3$ — > \/%tx becomes the SU(2) chiral Gross-Neveu S- 
matrix. This asymptotically free model has to correspond to one of the possible 
RG trajectories in the two-parameter SG language. It is easy to see that it has to 
be the Q = trajectory, since this is the only trajectory going through the origin 
(5q = «o = 0) of the parameter space. More precisely, the chiral Gross-Neveu model 
must correspond to the negative half of the Q = trajectory, which is an UV 
asymptotically free trajectory. Making the identification 

So = --g 2 , (17) 

7T 

where g is the coupling of the SU(2) Gross-Neveu model, the Gross-Neveu beta- 
function becomes 

(3(g) = --g 3 + ^g 5 + ... (18) 
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Using the results of Amit et al., (|), B = 2, using the results of |J and ||], (|7|), 
B = 4/3 and finally B = 8/3 if we trust ||, Comparing ((T8|) to the results of the 



beta-function calculations performed directly in the fermion language [BJ we see that 
the correct Gross-Neveu beta-function is reproduced if B — 2. Thus the two-loop 
results of Amit et al. are correct after all! This was the observation^ that served as 
our motivation for the present study. The correctness of the two-loop Gross-Neveu 
beta-function coefficient has been checked by studying the system in the presence 
of an external field 0]. Using this method the value of this coefficient can be read 
off from the bootstrap S-matrix and the results are in agreement with |J. 

We now turn to the explicit calculation of the renormalization parameters (|4]|J). 
The first quantity we consider is the two-point function of the U{1) current 

(U*)UV)) = J jffiS (^T " <W) e**-til(p). (19) 

The advantage of considering this physical quantity is that it is IR stable. Putting 
m = we find 

I(p) = I J 1 + 5 + ^ (inpa + K + i )+ ^(lnpa + Kf + . . \ , (20) 



7i { 32 V 2/ 16 

where K = — r'(l) — 1 — In 2. Since the current is conserved there is no operator 
renormalization required here and fl2"0| ) must become finite after the substitutions 
(|4],[|). From this requirement we get 

h = -&, 9i = -2, h = -&- (21) 

To determine the remaining two-loop coefficients g2 and we have to calculate 
Z a , the renormalization constant corresponding to olq. For this purpose we need a 
quantity with a perturbative series starting at 0{o.q). We have chosen the 2n-point 
correlation function 

X=(A(x 1 )...A(x 2n )), (22) 

where 

A{z) = fy™ e%«*\ (23) 

Although, in contrast to the Noether current, the operator fl23|) needs to be renor- 
malized, for large enough n the dimension of (EST) is so small that there is no operator 



2 We thank P. Forgacs who made this observation first and called our attention to it. 
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mixing and the operator renormalization constant can simply be determined from 
the correlation function 



Y = (A( Xl ) . . . A{x n )A*{ yi ) . . . A*(y n )}. 
A second order calculation gives 



(24) 



X n 2 2n 4 64n 3 



+ L 



o ? 



On 



where 



64n 2 128 



M 



2n 4 ' 
W 



n 



64n 3 
+ W 



+ 



(25) 



rii<i 




Yik<i 


Vk - yi\ 


IL,A; 


Xi - Vk 





(26) 



L = In Ma n and the dots stand for finite 0(a 2 ] ) terms as well as higher order terms. 
W(ijl) is defined by 

r 1 r 1 dz 

W(/i) = -l+ dzz^F{fi,fi;l;z)+ — [F(jjl, //; 1; z) - 1 - fi 2 z] , (27) 
Jo Jo z 

where F(oe, (3; 7; z) is the standard hypergeometric function. (p5|) can be made finite 
by the renormalization Yr = Z 2n Y, where 



n ^ 2n 2 



1 »v + ^V + M« 2 + ... 



with 



fci 



1 n 
+ 



64n 128 



64n 



W -)+W — 



n 



n 



(28) 



(29) 



For the 2n-point function X a second order calculation gives 



IQty 
+ 



not? 



128n + 64 



* + 4 H nW 

n 



(30) 



+ . . 



where 



7v=n 



i<j 



d 2 z- 



IL \ z - x i 



(31) 
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and 



* = -1 + -I In (Na- n( - 2n -^) - — Yf d 2 z ^IV'l 
^ 2 V > nF^J UAz-x^ 



(32) 



In (|30|) the dots represent finite terms of O(«o) and 0(5q) as well as higher terms. 
Renormalizing X by requiring Xr = Z2 n X to be finite after coupling constant 
renormalization gives 



23 = 



and 



92 



1 

16 



n 

128 



W 



n 



W 



(33) 



At first sight #2 seems to be n-dependent which would mean that the 2n-point 
function (^) cannot really be made finite with wave function plus coupling constant 
renormalization. Luckily, however, one can see that using the identity 



W{/i) - W(-fj) = -2/j, 



(H<i) 



(34) 



satisfied by the hypergeometric function, g2 is equal to the n-independent constant 
—5/64. Moreover, (|31f ) together with (PT[) reproduce (IS), the results of 0. The 
nontrivial cancellation of the n-dependence makes us more confident that these are 
the correct two-loop coefficients. 
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